利用AMMI模型分析寒地水稻3个品质性状的基因型与环境互作
刘丽华1, 胡远富2, 陈乔1, 李红宇1, 钱永德1, 吕艳东1, 郑桂萍1, 左豫虎1,*
1黑龙江八一农垦大学农学院, 黑龙江大庆163319
2黑龙江省农垦牡丹江管理局双峰农场, 黑龙江密山158308
* 通讯作者(Corresponding author): 左豫虎, E-mail:zuoyhu@163.com
摘要

为给寒地水稻品种审定和推广应用提供科学依据, 将黑龙江省垦区6个寒地水稻品种(系)种植在6个不同生态点, 利用AMMI模型对其营养食味品质性状(蛋白质含量、游离脂肪酸含量、直链淀粉含量) 进行了稳定性和适应性分析。结果表明, 蛋白质含量、游离脂肪酸含量、直链淀粉含量3项指标在基因型间、环境间及基因型×环境互作间的方差均达到极显著水平, 这3项指标的交互效应主成分值(IPCA)差异也达到显著水平; 影响营养食味品质稳定性的主要环境气候因子是穗期(抽穗至成熟期)日平均气温和日均气压, 穗期温度高、光照充足和少雨的环境有利于寒地稻米良好品质的形成; 6个参试品种(系) 营养食味品质稳定性为建07-1023>空育131>龙粳20>饶选06-06>北粳9005=垦稻08-924, 6个农场对水稻品质影响为查哈阳农场>梧桐河农场>853农场=军川农场>850农场>290农场。

关键词: 水稻; 寒地; AMMI模型; 营养食味品质; 稳定性; 适应性
Interaction of Genotypes with Environments for Three Quality Traits of Rice in Cold Region by AMMI Model
LIU Li-Hua1, HU Yuan-Fu2, CHEN Qiao1, LI Hong-Yu1, QIAN Yong-De1, LÜ Yan-Dong1, ZHENG Gui-Ping1, ZUO Yu-Hu1,*
1 Heilongjiang Bayi Agricultural University, Daqing 163319, China
2 Shuangfeng Farm, Mudanjiang Bureau of Heilongjiang Land-Reclamation, Mishan 158308, China
Abstract

In order to provide a scientific basis for varietal approval and application of rice in cold region, the genetic stability and adaptability of six rice cultivars from Heilongjiang Reclamation region planted in six different ecological sites in Heilongjiang Province were evaluated for nutritional and eating quality traits (protein, free fatty acid and amylose contents) based on AMMI (additive main effects and multiplicative interactions) model. The results showed that differences in nutritional and eating quality traits among genotypes, environments and among genotype × environment interactions were all significant at 1% probability level. The protein content, free fatty acids content and amylose content were also significantly different in IPCA value. Environmental factors affecting stability of nutritional and eating quality traits were daily mean air temperature and daily mean air pressure from heading to maturity stages. The environment with high temperature, enough illumination and little rain at rice heading stage was advantageous to produce good quality of rice in cold region. The order of stability of nutritional and eating quality traits for the six cultivars was as follows: Jian 07-1023>Kongyu 131>Longjing 20>Raoxuan 06-06>Beijing 9005=Kendao 08-924. The order for identifying rice quality of six test farms was Chahayang Farm >Wutonghe Farm >853 Farm =Junchuan Farm >850 Farm >290 Farm.

Keyword: Rice; Cold region; AMMI model; Nutritious and eating quality; Stability; Adaptability
引言

稻米品质性状多表现为易受环境影响的质量-数量性状[1,2], 这给优质水稻的育种和推广利用带来一定的影响。众多稻米品质性状中, 营养食味品质是稻米品质改良的最终目标。水稻营养食味品质的提高和优质品种的高效化利用都涉及生态环境和栽培技术措施等多项因素。因此, 研究不同环境条件下稻米品质的品种稳定性, 发掘对环境钝感的优质水稻品种(系), 是加快优质水稻品种培育与推广应用的有效途径。

在作物品种区域试验中, 不同品种在不同地点的产量和品质差异取决于品种的基因型和环境互作效应。要对品种全面合理地评价, 除了需要准确、可靠和有代表性的试验资料外, 还离不开合理有效的试验分析模型和方法。前人在研究作物品种稳定性和基因型与环境互作分析方面做了大量工作, 提出和分析了众多的稳定性分析统计模型和方法[3,4,5,6,7,8]。其中, AMMI模型把方差分析和主成分分析结合在一起, 充分利用试验所获得的信息, 最大程度地反映互作变异, 通过从加性模型的互作项中进一步分离出若干个乘积项之和, 提高估计的准确性, 相比其他方法, 该方法得出的结果更可靠, 已被广泛应用于不同作物品种区试评价[9,10,11,12,13,14,15,16], 但AMMI模型用于高寒地区水稻产量和品质的研究较少。本研究利用AMMI模型对2011年黑龙江省垦区寒地水稻区试品种营养食味品质进行分析, 以参数 Di评价参试品种的稳定性,参数 Dj估计试点对品种的影响, 分析品种和试点互作效应评价品种对试点的特殊适应性, 以期为寒地水稻品种审定和推广应用提供科学依据。

1 材料与方法
1.1 试验材料

应用2011年黑龙江省垦区寒地水稻区域试验数据。参试品种6个, 试点6个, 参试品种和试验地点在黑龙江省垦区均比较具有代表性。试验品种(系)和地点情况及代码见表1

1.2 试验方法

1.2.1 试验设计 采用随机区组设计, 3次重复, 小区面积为30 m2, 行距和穴距分别为30 cm和12cm, N、P、K施肥水平和田间管理均按黑龙江省农垦水稻区域试验方案统一实施。

1.2.2 测试项目与方法 将采集的试验样品于阴凉通风处风干至含水量14.5%, 用FC-2K型实验砻谷机(YAMAMOTO, 离心式)加工成糙米, 用FOSS近红外谷物品质分析仪(Infratec 1241)测定其蛋白质、脂肪和直链淀粉含量, 这3个性状可以代表主要的水稻营养食味品质性状。

1.2.3 统计方法 在基因型与环境互作效应显著的基础上, 进行同质性测验和联合方差分析, 利用AMMI模型进行品种稳定性分析。

[14]

式中, Yger代表第 g个基因型在第 e个环境中第 r次重复观测值, u是总体平均值, αg是基因型平均偏差(即各个基因型平均值减去总平均值), βe是环境的平均偏差(即各个环境的平均值减去总平均值), λn是第 n个主成分分析特征值, γgn是第 n个主成分基因型主成分得分, δen是第 n个主成分环境主成分得分, n代表模型主成分分析中主成分因子轴的总个数, θger代表误差。 为所估算的基因型与环境交互作用(G×E), λn0.5 rgn λn0.5 δen分别代表基因型G和环境E交互作用的第 n个交互作用主成分(IPCAn)。在所有显著的IPCA上有较小值的基因型或环境就为稳定的基因型或环境, 因此, 在IPCA双标图上越接近坐标原点的基因型或环境越稳定[3]

参照吴为人[17]的方法计算品种稳定性参数 Di。它是指一个品种(或基因型)在交互作用的主成分(IPCA)空间中的位置与原点的欧氏距离。

表1 试验品种(系)和试验点情况及代码 Table 1 Tested varieties(lines), test sites and code

式中, n是显著的IPCA个数, γin是第 i个基因型在第 n个IPCA上的得分, ωn是权重系数, 它表示每个IPCA所解释的平方和占全部IPCA所解释的平方和的比例。用 Di为所有基因型给出相应的定量指标, 品种的 Di值越小, 其稳定性越好。

采用Microsoft Excel 2000和DPS 7.05[18]分析 数据。

2 结果与分析
2.1 基因型与环境互作的AMMI 模型分析

表2可以看出, 本试验中各个参试品种的品质性状在不同试验区表现不同的平均值, 说明水稻品种间同一品质性状的稳定性不同。平均蛋白质含量以龙粳20最高, 为8.08%, 垦稻08-924其次, 为7.73%, 建07-1023最低, 为7.48%; 平均游离脂肪酸以垦稻08-924最高, 达21.47%, 空育131其次, 为21.41%, 饶选06-06最低(17.03%); 平均直链淀粉含量以垦稻08-924最高, 达11.02%, 饶选06-06最低, 为9.91%。

同一试验区不同品种的3个品质性状平均值不同(表3), 说明不同环境对水稻同一品质性状的影响不同。平均蛋白质含量以290农场最高, 为8.07%,军川农场其次, 为7.97%, 850农场最低, 为7.21%; 平均游离脂肪酸含量以查哈阳农场含量最高, 达21.49%, 军川农场其次, 为20.57%, 850农场最低, 达18.39%; 平均直链淀粉含量以853农场最高, 达11.57%, 梧桐河农场最低, 为9.31%。

方差同质性检验表明, 各试验点数据误差同质, 在此基础上对3个稻米品质指标进行联合方差分析(表4)。基因型、环境及其互作平方和占处理平方和的百分比(即SS%)大小可以反映它们对不同品质性状影响的大小。表4可以看出, 环境对蛋白质和直链淀粉的形成影响最大, SS%分别为53.56%和56.57%; 基因型对游离脂肪酸的形成影响最大, SS%达到52.75%; G×E对试验中3个品质指标影响均较大。从表5可以看出, 3个品质指标在品种(系)间、环境间差异及G×E互作效应均达到极显著水平, 故有必要利用AMMI模型对这3个营养食味品质性状进行稳定性分析。

2.2 营养食味品质稳定性及地点影响分析

AMMI模型稳定性分析表明(表4), 蛋白质含量和直链淀粉含量主成分因子IPCA1、IPCA2的F值达极显著水平; 游离脂肪酸含量IPCA1的F值达极显著水平, IPCA2 F值达显著水平。蛋白质含量、游离脂肪酸含量和直链淀粉含量的G×E第1和第2个交互作用的主成分之和(IPCA1+IPCA2)已分别解释G×E总变异平方和的84.68%、81.84%和76.74%。因此, IPCA1、IPCA2代表的互作部分, 能对营养食味品质稳定性做出判断。表5列出了6个品种(系)和6个地点的IPCA1、IPCA2及相应稳定性参数 Di值和Dj值。6个供试品种(系)中蛋白质含量稳定性最好的是建07-1023, 其次是饶选06-06, 稳定性最差的是龙粳20和北粳9005; 游离脂肪酸含量稳定性以龙粳20最优, 北粳9005最差; 北粳9005直链淀粉含量最稳定, 其次是建07-1023, 再次是饶选06-06。

表2 同一品种在不同试验区的3个品质性状平均值 Table 2 Average value of three quality traits of same varieties from different test sites (%)
表3 同一试验区不同品种的3个品质性状平均值 Table 3 Average value of three quality traits of different varieties from same test sites (%)
表4 稻米营养食味品质性状的基因型和环境互作效应分析 Table 4 Analysis of the interaction of genotype and environment for rice nutritional and eating quality traits

Dj值可以反映各试验地点对品种各营养食味品质的影响, Dj值越大, 说明影响越强, 反之, 越弱, 影响强的农场更适宜作为区试地点。由表5可以看出, 本试验6个试点中, 对蛋白质影响最强的是梧桐河农场, 其次是853农场, 影响最差的是290农场; 对游离脂肪酸影响最强的是850农场, 其次是查哈阳农场, 影响较差的是290农场; 对直链淀粉影响最强的是军川农场, 其次是查哈阳农场, 影响最差的是290农场。

2.3 品种适应性分析

由于蛋白质、游离脂肪酸和直链淀粉含量的G ×E前2个主成分之和已能解释70%以上的总变异(表4), 故以IPCA1为X轴, IPCA2为Y轴建立的AMMI双标图(图1)也能较好地反映品种的稳定性和地点的影响。在AMMI双标图上越接近坐标原点品种(系)稳定性越好, 越远离坐标原点的地点影响越强, 双标图反映的品种稳定性和地点的影响大小与表3的分析结果是一致的。

表5 参试品种与地点蛋白质、游离脂肪酸和直链淀粉含量的稳定性参数 Table 5 Stability parameter of protein, free fatty, acid and amylose content of the varieties from different test sites

图1 寒地水稻营养食味品质性状的AMMI交互作用双标图Fig. 1 AMMI biplots ( IPCA1- IPCA2 ) for nutritional and eating quality traits of rice in cold region

品种在地点图标与原点连线二维空间的垂直投影代表其在此试点的最大交互效应, 在正向连线上的最大投影代表此品种(系)在此试点表现出最佳适应性, 若垂直投影在地点和原点反向延长线上则表现不适应性。蛋白质项中, 龙粳20在查哈阳农场、空育131在梧桐河农场、北粳9005在853农场, 这些图标与原点的连线有较大的垂直投影, 表明这些品种(系)在这些试点表现出最佳适应性; 游离脂肪酸项中, 空育131在850农场、建07-1023和垦稻08-924在查哈阳农场、北粳9005在853农场均表现出最佳适应性; 直链淀粉项中, 龙粳20在850农场、饶选06-06在查哈阳农场、垦稻08-924在梧桐河农场均表现出最佳适应性。

上述分析表明, 建07-1023的蛋白质、龙粳20的游离脂肪酸、北粳9005的直链淀粉含量受环境影响最小, 适应性最好; 此外, 其他参试品种蛋白质、游离脂肪酸和直链淀粉含量在活动积温较高、水稻抽穗后日照时间长、降水量少的地点适应性较好。

2.4 G×E互作及其与环境气候因子的相关性分析

表6列出了试验点的各项IPCA值与穗期的环境气候因子的相关系数。蛋白质含量项分析显示, 其稳定参数 Di与日均气温达到显著正相关, 相关系数为0.82*, 与日均降雨量呈负相关, 但未达显著水平, 结合图1的分析结果, 龙粳20在查哈阳农场、空育131在梧桐河农场、北粳9005在853农场的G×E互作中, 穗期日均气温较高、日均降雨量较低的环境下有较好稳定性。同理, 对游离脂肪酸而言, 其稳定参数 Di与日均气压达到显著负相关, 相关系数为0.89*, 与其他环境气候因子相关均不显著, 空育131、建07-1023、垦稻08-924、北粳9005在日均气压较低、平均日照时数较多的地点均表现出最佳适应性。对直链淀粉而言, 稳定参数 Di与日均气温、日照时数呈负相关, 与日均降雨量和日均气压呈正相关, 但相关性均未达显著水平。

3 讨论
3.1 AMMI模型及品种稳定性分析

培育优质而稳定的品种是水稻品质育种最重要的目标, 在分析品种稳定性的众多数学模型中, AMMI模型最佳[7]。本试验使用AMMI模型分析发现, 环境对供试寒地水稻品种蛋白质和直链淀粉含量影响最大, 其次为G×E, 基因型最小; 对游离脂肪酸的影响为基因型>G×E>环境; 本试验中供试品种3项品质指标的综合稳定性呈建07-1023>空育131>龙粳20>饶选06-06>北粳9005=垦稻08-924的趋势, 试点影响为: 查哈阳农场>梧桐河农场>853农场=军川农场>850农场>290农场。说明寒地水稻品种及营养食味品质指标在不同基因型间、地点间及互作方面均存在差异, 这与张坚勇等[12]和万向元等[19]的研究结果一致。虽然稳定性参数 Di Dj可定量描述品种的稳定性以及地点的适应性, 但品种间稳定性差异以及地点间适应性差异是否达到显著水平, 还有待进一步研究。

表6 3个品质性状稳定参数与气候因子之的相关系数(抽穗至成熟期) Table 6 Correlation coefficients between stability parameter of three quality traits and climate factors from heading to maturity stage

依据AMMI分析的前两位主成分值及其相应的稳定性参数 Di值, 本试验对参试品种(系)营养食味品质与环境间的稳定性评价结果与现实种植中品种的稳定性表现基本一致。综合比较表明, 建07-1023营养食味品质最稳定; 空育131和龙粳20营养食味品质较稳定; 饶选06-06、北粳9005、垦稻08-924三个品系品质表现较差。现实生产中, 空育131和龙粳20两个品种的生育期基本一致, 品质稳定性均较好, 至今一直作为品比试验和区试试验中的对照品种; 各品系的营养食味品质还值得进一步探讨。

3.2 品种营养食味品质适应性分析

营养食味品质属数量性状, 除受遗传因素控制外, 还受环境条件的影响[20,21]。本研究结果显示, 6个寒地水稻品种(系)在不同地点种植, 营养食味品质适应性存在差异, 适应性因品种和环境而异, G×E之间存在很强的互作效应。建07-1023的蛋白质、龙粳20的游离脂肪酸、北粳9005的直链淀粉含量受环境影响最小, 适应性最好, 此外, 其他参试品种蛋白质、游离脂肪酸和直链淀粉含量在活动积温较高、水稻抽穗后日照时间长、降水量少的地点适应性较好。因此, 在品种研究中根据育种目标, 培育品质稳定性较好品种很重要, 但如果能发现品种和环境互作很强的品种也很有意义。AMMI模型是分析品种的多地点试验数据十分有效的工具, 但是目前只能对单个性状进行分析, 如何对多个性状进行综合分析还有待于进一步深入研究。

3.3 G×E互作及其与气候因子的相关性分析

G×E互作与气候因子相关分析表明, 不同试验抽穗期到成熟期的环境因子对蛋白质、游离脂肪酸和直链淀粉含量有较大影响(表4)。对蛋白质含量影响较大的气候因子是日平均气温, 两者呈显著正相关, 蛋白质含量与日均降水呈负相关; 对游离脂肪酸影响较大的是日均气压, 两者呈显著负相关, 游离脂肪酸与日均温度和日照时数均呈正相关; 直链淀粉与日均气温、日照时数呈负相关, 与日均降雨量和日均气压呈正相关, 说明穗期温度高、光照充足和少雨的环境有利于稻米品质的提高。已有研究指出, 一般品种的直链淀粉含量主要受灌浆期温度影响, 较低的生长温度能诱导粳稻的Wx基因的表达和GBSS酶活性的增加, 从而使直链淀粉含量增加, 口感变差[19], 这与本试验“直链淀粉稳定性与日均气温呈负相关”结果一致, 但相关性未达显著水平, 这可能是试验误差造成的, 有待进一步研究。

4 结论

基因型、环境及其互作对寒地水稻蛋白质含量、游离脂肪酸含量、直链淀粉含量影响不同, 其中G×E的影响较大; AMMI模型可以判断6个参试品种(系)营养食味品质的稳定性和6个试点的适应性, 本试验条件下, 建07-1023品质稳定性最好, 查哈阳农场最适合作物品质鉴定; 穗期温度高、光照充足和少雨的环境有利于寒地稻米良好品质的提高。

The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。

参考文献
[1] Zhang J-Y(张坚勇), Wan X-Y(万向元), Xiao Y-H(肖应辉). . Analysis on stability of eating quality of cooked rice (Oryza sativa L. ). Sci Agric Sin(中国农业科学), 2004, 37(6): 788-794 (in Chinese with English abstract) [本文引用:1]
[2] Miura K, Lin S Y, Yano M, Nagamine T. . Mapping quantitative trait loci controlling seed longevity in rice (Oryza sativa L. ). Theor Appl Genet, 2002, 104: 981-986 [本文引用:1] [JCR: 3.297]
[3] Piepho H P. . Analyzing genotype-environment data by mixed models with multiplicative effects. Biometrics, 1997, 53: 761-766 [本文引用:2] [JCR: 1.827]
[4] Wu Y-Q(吴元奇), Pan G-T(潘光堂), Rong T-Z(荣廷昭). . Study progress in crop stability. J Sichuan Agric Univ(四川农业大学学报), 2005, 23(4): 482-489 (in Chinese with English abstract) [本文引用:1] [CJCR: 0.5871]
[5] Hu B-M(胡秉民), Geng X(耿旭). Crop Stability Analysis (作物稳定性分析). Beijing: Science Press, 1993. pp68-70(in Chinese) [本文引用:1]
[6] Piepho H P. . Stability analysis using the SAS system. Agron J, 1999, 91: 154-160 [本文引用:1] [JCR: 1.794]
[7] Zhang Q-Y(张群远), Kong F-L(孔繁玲). . Comparison of statistical models for regional crop trial analysis. Sci Agric Sin(中国农业科学), 2002, 35(4): 365-371 (in Chinese with English abstract) [本文引用:2]
[8] Hu X-Y(胡希远), You H-L(尤海磊), Song X-F(宋喜芳). . Comparison of different models for crop stability analysis. J Triticeae Crops(麦类作物学报), 2009, 29(1): 110-117 (in Chinese with English abstract) [本文引用:1] [CJCR: 0.8529]
[9] Liu W-J(刘文江), Li H-J(李浩杰), Wang X-D(汪旭东). . Stability analysis for elementary characters of hybrid rice by AMMI model. Acta Agron Sin(作物学报), 2002, 28(4): 569-573 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.8267]
[10] Shi W-X(施万喜). . Analysis on high yield and stability of winter wheat cultivars in dryland of east Gansu by AMMI model. Agric Res Arid Areas(干旱地区农业研究), 2009, 27(3): 37-43 (in Chinese with English abstract) [本文引用:1]
[11] Ji J-H(冀建华), Liu G-R(刘光荣), Li Z-Z(李祖章), Liu Y-R(刘益仁). . Effects of long-term fertilization on stability of double cropping rice yield based on AMMI model. Sci Agric Sin(中国农业科学), 2012, 45(4): 685-696 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.4522]
[12] Zhang J-Y(张坚勇), Xiao Y-H(肖应辉), Wan X-Y(万向元). . Stability analysis for appearance qualities of rice cultivar. Acta Agron Sin(作物学报), 2004, 30(6): 548-554 (in Chinese with English abstract) [本文引用:2] [CJCR: 1.8267]
[13] Jiang K-F(蒋开锋), Zheng J-K(郑家奎), Zhao G-L(赵甘霖). . Stability of grain yield traits and their correlation in hybrid rice. Chin J Rice Sci(中国水稻科学), 2001, 15(1): 67-69 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.4083]
[14] Su Z-X(苏振喜), Zhao G-Z(赵国珍), Liao X-H(廖新华). . Stabi- lity of eating quality traits of japonica soft rice from different locations at different altitudes in Yunnan Province, China. Chin J Rice Sci(中国水稻科学), 2010, 24(3): 320-324 (in Chinese with English abstract) [本文引用:2] [CJCR: 1.4083]
[15] Yao X(姚霞), Li W(李伟), Yan Z-H(颜泽洪). . Application of AMMI model on the yield trait analysis of wheat regional trials. J Triticeae Crops(麦类作物学报), 2005, 25(6): 103-107 (in Chinese with English abstract) [本文引用:1] [CJCR: 0.8529]
[16] Jiang K-F(蒋开锋), Zheng J-K(郑家奎), Zhao G-L(赵甘霖). . Analysis of combining ability based on AMMI model. Acta Agron Sin(作物学报), 2000, 26(6): 959-962 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.8267]
[17] Wu W-R(吴为人). . An improvement on the method of variety stability analysis based on the AMMI model. Hereditas(遗传), 2000, 22(1): 31-32 (in Chinese with English abstract) [本文引用:1] [CJCR: 0.9787]
[18] Tang Q-Y(唐启义), Feng M-G(冯明光). Practical Application of Statistics Analysis and DPS Data Processing System (实用统计分析及其DPS数据处理系统). Beijing: Science Press, 2002 (in Chinese) [本文引用:1]
[19] Wan X-Y(万向元), Hu P-S(胡培松), Wang H-L(王海莲). . Analysis on stability of AC, GT and PC in rice varieties (Oryza sativa L. ). Chin J Rice Sci(中国农业科学), 2005, 38(1): 1-6 (in Chinese with English abstract) [本文引用:2] [CJCR: 1.4083]
[20] Shen X-H(沈希宏), Yang S-H(杨仕华), Xie F-X(谢芙贤). . Genotype × environment interaction and its relationship to climate factors in rice regional trial. Chin J Rice Sci(中国水稻科学), 2000, 14(1): 31-36 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.4083]
[21] Liu J(刘健). . Research advances on the effects of environmental factors on rice grain quality. J Hubei Agric Coll(湖北农学院学报), 2002, 22(6): 550-553 (in Chinese with English abstract) [本文引用:1]